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I. INTRODUCTION. 



The main aim of this paper is to demonstrate the method caUed "the Bosonization of Nonlocal 
Currents" (BNC), used for calculations of bound states in a quark model, within the simplest 
relativistic quantum field model of two scalar fields with the Yukawa type interaction: 

L = $+(□ - M^)<E> + ^(f>{D - m^)(j) + g(f)^+^. 

A second aim is to clarify the relation between BNC and two widely used methods, employed in 
■ recent particle physics to calculate bound states of interacting particles, based on the nonrelativistic 

' Schrodinger equation (the S-method), and the relativistic Bethe-Salpeter equation (the BS-method), 

^\ , and to determine the conditions on parameters of a quantum field model dictating a definite method 

to be applied. 

It is shown that all these methods can be applied only in the weak coupling regime (the effective 
dimensionless coupling constant A should be less than 1). The basic parameter separating the 
relativistic and nonrelativistic pictures is ^ = namely, <^ 1 with A <c: 1 leads to the potential 
, picture, i.e., the bound state is described by the nonrelativistic Schrodinger equation. For ^ > 1 and 

04 ■ A < 1 the Schrodinger potential picture is not valid and the Bethe-Salpeter equation or the BNC 

should be used, where the BNC method has a slightly bit wider region of applicability. 
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On ' Great efforts fiave been made to understand liow bound states arise in the formalism of quantum field theory and to 
\ work out effective methods to calculate all characteristics of these bound states, especially their masses. Unfortunately, 
• we observe that there is no well-defined unique method, like the Schrodinger equation in nonrelativistic quantum 
rS i' mechanics, which can be used practically for any problems of nonrelativistic quantum physics. We can conclude that 
JL QFT of today is not well suited to describe bound state problem, (see, for example, 

The analysis of a bound state is simplest when the constituent particles can be considered to be nonrelativistic, 
(~| , i.e., when they travel at speeds considerably less than c. The physical evident criterion to tell that a bound state is 
ILJ • nonrelativistic is for the binding energy to be small compared to the rest energies of the constituents. The theoretical 
criterion is that the coupling constants should to be weak and masses of intermediate particles (photon in QED, 
mesons in nuclear physics, gluons in QCD) should be small in comparison with masses of constituents. 
5-H ' The best example is hydrogen-like systems, which can be considered nonrelativistic, and the experiments with great 
accuracy were required to develop the methods to calculate next relativistic corrections (see, for example, fll-^). 

The situation in nuclear and particle physics is completely different. First, the coupling constants are not small 
any more. Second, in nuclear physics, although the binding energy is relatively small in comparison with nucleon 
masses, the masses of intermediate mesons realizing the strong nuclear interaction are not small. In particular, the 
most adequate description of the deuteron can be done by the Bethe-Salpeter equation where the contribution of all 
light mesons should be taken into account (see, for example, Q). 

In particle physics, only hadrons made out of heavy quarks can be considered by the nonrelativistic potential meth- 
ods although the decays of heavy hadrons into light ones require relativistic methods to describe these transformations. 
The most familiar light-quark states are intrinsically relativistic, so that they require pure relativistic methods. Be- 
sides, they are constituted at distances where the confinement phenomenon should be taken into account. In addition, 
one may ask whether the free Dirac equation applicable to describe the light quarks in this region? Therefore, using 
the nonrelativistic Schrodinger equation to describe the light-quark systems by fitting parameters of potentials (see, 
for example, [^|j^]) can be considered heroic attempt to understand light meson physics by unsuitable methods in a 
very rough theoretical approximation. 

In relativistic quantum field theory, bound states are identified by the occurrence of poles of corresponding ampli- 
tudes or Green functions with appropriate quantum numbers. These poles have a nonperturbative character, so that 
they can arise as a result of a nonperturbative rearrangement of series over a coupling constant. The investigation of 
nonperturbative properties was done by establishing integral equations among amplitudes and Green functions, using 
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the specific structure of a Lagrangian. One should say that these equations, having absolutely general form, in reality 
can be used when the kernels contain contributions of the lowest Feynman diagrams only. It implies that in some sense 
the coupling constant should be small enough. The Bethe-Salpeter equation is the most important integral equation 
of this type and it is widely used, especially for calculation of relativistic corrections in hydrogen-like systems and 
deuteron physics (see, for example, [0J|,0,^)- 

Now we would like to pay attention to the so-called Z2 = approach (see, for example, [p|jic[|) which is not used so 
widely, although the bosonization of QCD introducing bilocal boson-type fields was developed in The objection 
and prejudice against this method are based on the persuasion that in local quantum field theory we should have local 
interaction only. For example, the pion as a quark-antiquark bound state is represented by the local term 77(575^). As 
a result the renormalization constant Z2 for the pion contains the ultraviolet divergence and therefore the condition 
Z2 ~ makes no practical sense. The answer is that the local vertex does not lead to any bound state. The vertex 
connected with the formation of a bound state should be nonlocal. Moreover, this nonlocal vertex in the nonrelativistic 
limit is directly connected with the nonrelativistic wave function of this bound state. The question arises how to find 
this nonlocal vertex. In the approach ||ll[] this vertex is a solution of the Bethe-Salpeter equation, which can only be 
solved by difficult numerical computations. 

Recently, in papers jl^ we began to develop the Model of Induced Quark Currents based on the assumpsion that 
(anti-)self-dual homogeneous gluon field realizes the physical QCD vacuum. To calculate the mass spectrum of 
bosons as quark-antiquark bound states, the method of so-called Bosonization of Nonlocal Currents (BNC) was used. 
Numerical results for the boson spectrum turn out to be in good agreement with the experimental ones. This method 
is quite close to the Z2 = method and to the bosonization of QCD introduced by A similar approach was 

developed for fermionic many-body problems in statistical physics in [^. The idea of BNC consists of two point. 
First, we write the Green function containing a bound state in the functional integral representation in which the 
Gaussian measure being the lowest order of the bosonization of bilocal currents contains the Bethe-Salpeter kernel 
written in symmetric hermitian form. Second, we have guessed the analytical form of the orthonormal system of 
functions very close to the eigenfunctions of the Bethe-Salpeter equation. 

Now we would like to attract attention to the following quite important point which is not usually stressed clearly 
enough. Because exact solutions for any realistic quantum field model are not known one or another approximation 
should be used. In any case, we should have possibilities, at least in principle, to evaluate the theoretical accuracy 
of a chosen approximation. The most simple, although maybe not the most accurate, way is to require the effective 
coupling constant to be noticeably smaller than unity. The perturbation series in the coupling constant is the most 
popular method of calculations, especially in quantum electrodynamics or the theory of weak interactions. However, 
if a perturbation series is rearranged by a subsummation, the problem of the effective coupling constant becomes 
quite topical. Thus, if the theory is reformulated somehow, the effective coupling constant should be defined for each 
particular case and should be smaller than one. Only in this case we can trust our computations. 

The aim of this paper is to formulate the BNC-method and to clarify the connection between the principal ap- 
proaches used in atomic, nuclear and particle physics of today to calculate bound states of interacting particles, 
namely, the nonrelativistic potential Schrodinger equation (S-method), the relativistic Bethe-Salpeter equation in the 
one-boson-exchange approximation (the BS-method) and our BNC-method in the one-loop approximation. We answer 
this question by considering a relatively simple quantum field model. The results will provide a deeper understanding 
of approximations used in the well-known standard approaches. 

An example of the above-mentioned simple quantum field model is the Yukawa interaction of charged scalar bosons 
described by the field $ and neutral bosons described by the field (p. The Lagrangian density is 

L{x) = $+(□ - M2)$ + i0(n - rn^)(t) + 5$+$^. (1) 

This model is frequently used as the simplest pattern of QFT in many discussions, although this system is not stable 
from a strict point of view because the Hamiltonian is not bounded ^from below. This model has been investigated 
by various methods (see, for example, [p]J^Jl5[| and the recent paper and references therein). 

In this model it is possible to retrace all details of bound states arising in quantum field theory. Generalization 
to the case of the Dirac field presents no difficulties of principle and leads to technical problems connected with the 
algebra of 7-matrices only. This model is superrenormalizable so that the renormalization procedure has the simplest 
form. The main aim of this paper is to understand the general mechanism of bound states arising in this quantum 
field model and to clarify conditions on parameters of this quantum field model dictating a definite method mentioned 
above to be applied. 

The model contains three dimensionless parameters: 



where Mf, = 2M — AM is the mass of a supposed bound state, AM is the mass excess or the binding energy. 

Three parameters, A, ^ and 6, are not independent. The standard formulation of the problem is to find 6 if ^ and 
A are given. We formulate the problem in another way: what is the region of changing b for a fixed ^, if the effective 
coupling constant (which can differ ifiom A) is smaller than unity? 

The parameter ^ is supposed to be smaller than 1. Our aim is to find the condition under which the mass M(, of a 
bound state lies in the interval 

< Mb <2M or < 6 = (^^^ < 1, 

i.e. this bound state should be stable. 

It turns out that all approaches mentioned above require the effective coupling constant to be small enough. The 
value of the parameter ^ = plays the crucial role separating nonrelativistic and relativistic approaches. Namely 

• the Potential nonrelativistic picture takes place for 



and, therefore, 



1 - & 



\2M 



«1, 



or AM < 2M, 



i.e., the binding energy is very small; 
• the Bethe-Salpeter approach and Bosonization of Nonlocal Currents take place for 



m 

^=M^1' 



A = 



167rM2 



< 1, 



(4) 



and the binding energy can be up to AM < (.2 ~ .3) • 2M. 

The result of this paper is that all methods under consideration can be used in the weak coupling regime only, 
but the nonrelativistic potential picture takes place provided the mass of an exchange particle m to be very small in 
comparison with that of the constituent particle M. The BS-approach and the Bosonization of Nonlocal Currents 
can be used for any relations between the masses of exchange and constituent particles although the latter method 
has a somewhat wider range of applicability. 

From my point of view, it is just the BNC-formalism to be the most attractive method to study bound states in 
QFT. It provides to represent results in an analytic form and to evaluate the theoretical accuracy of approximations. 

This work has been supported by the Russian Foundation for Fundamental Research N 96-02-17435a. 



II. THE FORMULATION OF THE PROBLEM 



We consider the interaction of scalar charged scalar fields $(a;) and a neutral scalar field 4>{x). All the consideration 
is given in the Euclidean metrics. The total Lagrangian is done by (|l|). The propagator of a scalar particle with the 
mass M is 

f dp e^P'-^^y^ 

^M(--2/)=y (^-^^^^^ (5) 

The object of our interest is the four-point Green function 

G{xi,x2,yi,y2) (6) 
= J J D^D^+Dcp-^+{xi)^x2)<i>+{yiMy2)-e^^''^'''-'^\ 

£'[$+,$, 0] = / dxL{x), 
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where an appropriate normalization should be introduced. 

The four-point Green function (||) contains all information about possible bound states in the channels and 
The particles $ can be called constituent particles. We take an interest in bound states in the channel <!>"'"$. 
The quantum numbers Q can be fixed by an appropriate vertex 



$+(xi)$(x2) ^ Jq{x) = ($+Vq$), = $+(x)T/Q(p,)$(a;), 
dx — dx ■ 



1 

Px= - 



The nonlocal vertex Vq{Px) defines the quantum numbers Q of the current Jq = {^~^Vq^). This vertex can be 
represented like 



VQiPx) = J duVQ{u)e^-^-, (7) 
The current Jq — ($+Vq$) can be written as 

Jq{x) = ($+Fq<J>), = {<^+{x)VQ(Pxmx)) (8) 
du $+(a; + u)Vq{u)(^{x - u). 

The Green function with quantum numbers Q is defined as 

Gq{x ~y)^ j D<i>D^+ J JQ{x)JQ{y) e^t*^'*'^!. (9) 

If in this channel a stable bound state with the mass Mq = Mt < 2M does exist, the Green function Gq{x) has the 
following asymptotic behaviour 

Gq{x) - e-^'^«l^l for |x| ^ oo, 
so that the mass of the state Jq = ($+Vq$) can be found as 

Mb = Mq = - lim ^ In Gq {x) . (10) 

x^oo \x\ 

The problem is to calculate the functional integral in representation and find the mass M{, according to ( |l0[ ) . 

If we consider perturbation expansion over the coupling constant g for the four point Green function G{xi,X2,yi, 2/2), 
we will get a series of the Feynman diagrams describing an interaction of two particles $. This series can be written 
in the form of the Bethe-Salpeter equation (see, for example, Bound states of two particles <I> in a channel 

Jq{x) = {^^Vq^)x can be found as solutions of this equation. 

We proceed in another way. First, our aim is to obtain for the Green function (^ the functional integral repre- 
sentation in which the term being responsible for bound state creation would be written in the explicit form and for 
remaining corrections the smallness criterion would be defined. Second, we want to formulate the method of analytical 
calculations of binding energies and evaluate their theoretical accuracy. 

Fortunately, it is possible in the representation (^ to do the first integration either over the field ^{x) or field 
4>{x). Thus, we get two representations which arc the starting points of two approaches: the Potential picture and 
the Bosonization of Nonlocal Currents. 

I. Potential picture 

The integration in (||) over the charged scalar field <I> gives for the Green function (^: 

Gq ix-y)^ Gg') {x~y) + G'^f (x-y). (11) 

Here 
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G^q\x -v) = j du j dv j Dcf) 

■ \yQ{u)'DM{x + u,y + v|0)Vq(u)Dji/(2/ -v,x- , 
G^Q^x-y) = J duj dv J D(t) e^^I-^l 

• |VQ(u)DAf(.T + u,x-M|(/))| • ^VQ{v)'DM{y + v,y - v\(j))^ , 

where 
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Sp[4>\ = - dx (t){x){a - m^)(/)(x) - trln[l - g(j)D 
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The Green function DM(a;, ?/|(/') satisfies the equation 

[-□ + M2 - g<j>{x)]nM{x, y\cf>) = Six - v) (12) 

with 

^l[{x,y\(l)) = DM(a;,y|(/)) = \:»M{y,x\(i)). 

Here the functions G'q\x) and g'q'\x) are said to be "potential" and "annihilation" Green functions, respectively. 
The approach based on the representation ( pi] ) will be called the Potential picture. 

II. Bosonization of Nonlocal Currents 

The integration in (^ over the scalar field cj){x) gives 

Gq{x -y)^ j j D^D'^+ ■ e^^-[*+.*l . {^+VQ^),{^+VQ^)y, (13) 
where an appropriate normalization is implied and 

($+1?^/$) ^ j dx $+(x)(-n + A-/2)$(a;), 

^($+$Z?,„$+$) = Jdxjdy ^+ixMx)D^{x - y)^+{yMy). 
The approach based on the representation ( p^ ) will be called t/ie Bosonization of Nonlocal Currents. 

III. THE POTENTIAL PICTURE. 

The starting point of the Potential picture is the representation (pi]). 

A. The Green function T)M{^,y\(l>) 

First, the charged loops should be neglected, so that the Green functions Gp{x — y) and Gq{x — y) are represented 
by 

G''^\x-y)^ J dauv[4>] '^Aiix + u,y + v\<j)) -BMix - u,y - v\(j)), (14) 

G^'^^x-y) =^ J dauv[4'] T)m{x + u,x - u\(f>) ■ T>M{y + v,y - v\(j)) 
dcr„„[(/)] = duV{u) dvV{v) D4) e-5('^'(^)cr„V)_ 



We would like to stress now that the neglect of loops of scalar particles $ presupposes the dimensionless coupling 
constant A to be small enough. 

The solution of (|l|) can be represented by the functional integral (see, for example, [p^ ): 

oo ^ a . ct \ 

l ) 

with the boundary conditions ^(0) — y, ^(q) = x and the normalization 







B. The Green function G''^\x) 

The function G^-^^ (x) after integration over ip has the form 

Gpix) = (^)'/ d^idE^ e^-+2W-+'^-, 
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I dE,{*} = I d^.VK) I da, e-^''' I D^.exp i-Jdr^-^ \ {*}, 

I J 

U = 1,2; ui = u, U2 = -y), 

6(0) = "!^, ^i(ai) = x + u, ^2(0) = -w, ^2{a2)^x-u. 

where for simplicity we put y — 0. 

Our task is to get the asymptotic behaviour of the fmictions G'^^\x) for asymptotically large x — — > 00. To 
this end, let us introduce the following variables: 

X (3j 

a 



^i(/3) = n/3 + 77i(/3), 6(/3) = n/3 + %(/3), n^ = ^. 

X 



Then, one can obtain 

I duV{u) I dvV{v) (16) 



M ' ^ 



Stt^x, 

00 00 
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where 



DriiDrj2 exp 



(17) 



dl3 



'7i(0) = w, rii{x)=u, 



m{^) = mix) = -u. 



The "potential interaction" is described by the two-point nonlocal functional 

W^^[?7i,?72] = Wii + 2Wi2 + W22. 



(18) 



df3idf32 An (n(/3i - P2) + m(Pi) - ■ 







The functional integral for Jx{si, S2', x) looks like the Feynman path integral in the nonrelativistic statistic quantum 
mechanics for the four dimensional motion of particles ?7i(/3) and ?72(/3) with "masses" Msi and Ms2 where /3 plays 
the role of the imaginary time or temperature. The interaction of these particles is defined by the nonlocal functional 
Wii + 2W\2 + W22 which contains potential W\2 and nonpotential W\\ + W22 interactions. It should be noted that 
the structure of the functional integral (|T^) reminds the polaron problem (see ||l7|] ). 
The asymptotic form of the function J^)(si, S2; x) looks like 

j(^)(si,S2;x) - exp{-a;£;(si,s2)}, (19) 

where E{si, S2) is the energy of the lowest bound state. The asymptotic behaviour of the functional G^^^x) as a; cxd 
is determined by the saddle point of the integrals over si and S2 in the representation (p^. Substituting expression 
( p^ ) into (|l^), one can get 



Mb = - lim -llnG(-P)(a;), 

x^oa \x\ 



(20) 



= mm 

(si,S2) 



M 

Y 



1 

Si 



Sl 



1 

S2 



S2 + E{si,S2) 



M [- + s]+ E{s, s) 



The main problem is to compute the functional integral (0). This computation can be done by the variational 
methods (see [|l^). We plan to calculate this functional integral applying the Gaussian equivalent representation 
method which was sucessfuUy used sucessfuUy for the polaron problem (see [|6|). 

Besides one can see that the representations ( p^ ) and ( pT| ) does not really feel the explicit form of the vertex Vq (u) 
although it should extract a bound state with definite quantum numbers Q. It means, in fact, that in the general 
case, i.e., for any value of the coupling constant g, the modern analytical methods, applied to the functional integral 
( |l6| ) , allow to calculate with reasonable accuracy the energy of the lowest bound state only. 

C. The Green function G'-'^^x - y) 
The function G^^^ (x) after integration over (/) has the form 



dT2D,n{^j{n)-^,{T2)), (j = l,2). 
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ai a.2 
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Here ^i(O) = -u, ^^{ai) ^ u, 6(0) = -w, 6(a2) = 
One can see that in the Umit x oo W12 and 



" M 




_8tt^x J 





so that no bound state arises in this case. Thus, the annihilation channel does not contain any bound states. In other 
words, intermediate pure boson states of particles 4> cannot produce any bound state. 



D. The Nonrelativistic Limit 



In this section we obtain the nonrelativistic limit c ^ 00 for the loop function Gp{x) in (|16|). Our task is to 
introduce the parameter c in an explicit form into Gp{x) and then go to the limit c ^ 00. To this end, let us restore 
the parameter c in our formulas. We have 



(21) 



X^ = (X4, x) — > (ct, x), X = = V C^t^ -I- x2 ct, 



(l,x/ct) ^ (1,0). 



The propagator D„i{x) becomes 



dk f dv 



(2n) 



oi(t)t+kx) 



m 



dk r dv e'(^*+'''') 
(27r)3 J 2^ ■ k2 + k2 

— 00 

1 e^'"' 



V 

1 - ^ • 



1 



c2 k2 + 



1 1 e^*"" 



47r 7' c2 ^ ' Stt K \ c"* 

The parameter -i defines the radius of the nonrelativistic Yukawa potential, therefore, we have to keep 

K — mc — const 

to be finite in the limit c — > 00. 

Let us come back to the Green function (p^. According to (|2^) it reads 



Gp{t) 



M 

47rf 



+ S2 



Jp{si, S2;t) 



(22) 



In the functional integral (|17|) let us introduce the new varibles (3j = ctj . It is convinient to represent the four- vectors 
r]j in the one and three component form 



r?,(/3) = (a,(T),b,(T)), (j = l,2) 



and 



ai(T) = $(r) + 



a2(r) = $(t) - 
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bi(r) = R(r) 



r(r) 



b2(r) - R(r) - 



r(r) 



with the boundary conditions 



R(0) = R(i) = 0, r(0) = 2u, r(t) = 2v 



In the hmit c ^ oo the integral for Gp{t) over si and S2 can be calculated by the saddle-point method. The saddle 
points are si = S2 = 1 and we get 



(23) 



where the functional integral looks 



n n -f- / dr[d?+di+b?+bi]+Wil+2Wi2 + lV22 

Jp{t) = / / DaiDa2 / / DhiDh2 e o 

t 

D^D(j) I I DRDr e « 

The functionals Wu and W22 describe the mass renormalization in the nonrelativistic limit and we omit them. 
According to (|2|) and (p2|) the functional W12 (p^, where the argument of _D-function is 



n{f3i - f32) + miP) ~ miP) = (c(ti - T2) + ai(Ti) - a2(T2), bi(Ti) - b2(T2)), 
acquires the form in the limit c 00 



2Wl2 = 



167rM- 



■// 



dri dT2 



5{ti - T2) + -5'{ti - T2)(ai(Ti) - a2(T2)) 



1 e""*" 



= X dT 







8c2 



(rr)^ 



(1 + Kr) 



r c 7' 



O { -0R, ^R2 

c 



where r = r(T) = |bi(r) — b2(r)|. 

After integration over <i>, and R in the functional integral for Jp one can get 



r(t)=2v 



Jp{t) 



Dr ■ exp ■ 



dr 



^r2(r) + C/(r(r)) 



r(0)=2u 



[/(r)) = -A- 



A 



r 2Af2c2 



Pi 



2 „-2Kr 



A^ e 



8Mcc2 r2 



(24) 



(25) 



where Mc = ^ and pj = Mcrj ■ 

Here we keep the main terms contributing to the potential and neglect terms describing and the nonlocal interaction 
and terms of the order 0(l/c^). 

One can see that this representation for Jp (t) coincides with the Feynman path integral in quantum machanics for 
the Green function (v, t; u, 0) 
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Mr: 



(26) 



r(t)=2v / t 

K{w,t]u,Q)^ j Drcxpl-Jdr 

r(0)=2u I 

Q 

Here '0Q(r) and Eq are eigcnfunctions and eigenvalues for the quantun number Q connected with the space R'' of 
the Schrodinger equation 



2Mr.. 



C/(r) 



Tpqif) ^ EqiPq{y) 



(27) 



where U{y) < is the attractive potential. 

As a result, the Green function Gp{t) for t — > oo behaves like 



Gp{t) = / ^" ^(u)^q(u) 



If we choose 



then 



and, finally, for large t we have 



V{n) = / dui V^(m4,u) =?/.Qj,(u), 



Gp(t) e 



where iSgg is the energy of the bound state of two nonrelativistic particles in the quantum state Qo arising due to 
the potential U{v). The mass of the bound state in the nonrelativistic approach is 



1 



Mc^ + Eq„+0 



(28) 



Thus, in the nonrelativitic limit the relativistic vertex V{u) is connected with the nonrelativistic eigenfunction 

Viu) ^ / du4 Viu4, u) = V{u) = V'q(u) (29) 



and the bound state mass of two scalar particles is a sum of their masses plus the binding energy Eq determined by 
the nonrelativistic potential interaction. 



E. The nonrelativistic Yukawa potential and relativistic corrections. 

Thus in the nonrelativistic limit, we have the Schrodinger equation with relativistic corrections: 

H = Hy + ^{Up + Us) + (4 



(30) 



where Hy is the Yukawa potential 



Mr: 



M 



2M„„ r ' ™' 2 



and Up and Us are the first and second lowest relativistic corrections in ( p4D 
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Our problem is to find restrictions on the parameters A and ^ = ^ for which the ralativistic corrections Up and Us 
can be neglected. We proceed in the simplest way. Wc use the variation function 

where ,s is the variational parameter, for the ground state. This approach gives us quite a good qualitative and even 
semi-quantitative estimation of the background energy. It is sufficient for our aim. The calculations give 

Ev ~ mm — ; — ; — — ^ = — — • mm {s — G ■ 



(^'1*) 4M s \ (l + s)2 



8Mc 3 + s ' 
where the parameter s is defined by the equation 

m s(3 + s) 

The bound state can exist if s > 1, i.e. 

G > 4, or A > — 

The contribution of relativistic corrections can be evaluated as follows: 

_ mUp\^) _ 1 /m\2 

^ (*|C/r|*) 4 ' \m) ' 1 + s' 

_ (^|;7,|^) _ A m {l + sf 

' (*|C/r|*) 8 ' M ' 2 + s ' 

Thus the nonrelativistic picture takes place if 

\Ey\ < 2M, 5p<^l, Ss<^l 
It is easy to see that these inequalities take place if 



M 

It is valid for s 2 3 and a (1 -j- 2) • m. 



TP 

^= — <1, A<1. (31) 



F. Relativistic incompleteness of quantum mechanics of two particles. 

Here we would like to pay attention to the Schrodinger equation which describes two nonrelativistic particles 

2mi 2m2 

where the potential is attractive. 

Let us pass to the center-of-mass system in a standard way 

^ mifi + TO2r2 ^ ^ ^ 

K= , r = ri— r2. 

mi + m2 

The Hamiltonian takes the form 

2 2 

2M 2/i ^ ' 

17111712 

M = mi+ m2, IJ: ■ 



mi + m2 
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The solution of the Schrodinger equation H"^ = E"^ can be written as 

where p is the momentum of the total system; ^/'(r) is an eigenfunction of the equation 



4'{r) = —eipir), 



and — e (e > 0) is an eigenvalue of a bound state. Then, the eigenvalue or the energy of the state 'i'{R, r) for p 7^ is 
From the physical point of view this energy has no reasonable meaning. Indeed, we should get 



2M 2Mphys 
Mphys = mi + m2 - A, ^ ^ ^ 

i.e., the interaction between two particles should give the mass excess. 

On the other hand, the latter formula can be obtained from the relativistic energy in the nonrelativistic limit 



2 



= (mi + m2 - AM)c2 + + 0{p^) 

[nil + m.2 — AM ) 

9 .9 f AM\ 



- (mi + m^y + — ^- r - AMc' + 0[—-]+ 0{p' 

2(mi + m2) \ M J 

and the mass of the bound state equals 



Thus, the nonrelativistic Schrodinger equation describing two nonrelativistic particles can be cosidered as a fragment 
of an relativistic equation describing the relativistic interaction of two particles. The "true" relativistic equation 
should contain terms describing both the motion of the center of masses and their relative motion. 

IV. BOSONIZATION OF NONLOCAL CURRENTS. 

The starting point of Bosonization of charged currents is the representation ([T^ ) . 

A. Bilocal currents. 

Let us consider the four-field term (^3|) and introduce the bilocal currents: 

^($+$i?„,$+$) 

J J ^^^"^^^ $+(2/i)*(2/i)£'m(2/i -2;2)*+(y2)«'(2/2) 

= Y J J ^y^^y^ J^iyi,y2)J{yi,y2) = y(-^^'^) 

where 
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Jiyuy2) = - ?y2)($+(yi)3'(y2)), 

J^{yi,y2) = J{y2,yi)- 



The next step is to use the Gaussian representation 



where A+{yi,y2) = A{y2,yi) and 

(A+A) = 
(A+ J) = 



DA e"2( 



dyidyz ^^(yi, y2)^(yi, ^2), 
dyidy2 A+{yi,y2)J{yi,y2) = (^-'^^) 



dyidya ya) V-Dm(yi - ^2)^(2/2) 

= ($+(^+/d:;:)$). 

We can always represent the product $+(xi)$(a;2)<i>+(a;3)<i>(x4) as a combination of currents 

g^^+{xi)^{x2)<i>+ix3)<i>{xi) g^J{xi,X2)Jix3, x^). 

We have 

G{X1,X2]X3,X4) 

„2j(^ ^_^7/ \ „-(<I>+-D-/*)-i(A+A)-g (A+J) 



D^D<i>+ J DA g\J{xi,X2)J{x3,X4) ■ e- 
Now we integrate by parts over A and and calculate the Gaussian integral over $ and <i>+: 

SIA] 



G{xi,X2;x3,Xi) = I I DA A{xi,x2)A{x3,Xi) e 
1 



S[A] = -^(A+A) - trln [l + giA+^/D^)DM 



Matrix operations are defined by formulas 

(A+y/D^) = A+{xi,X2)y/Dra{xi ~X2), 



iA+^/D^)D 



M 



= / dy A'^{xi,y)\jD„i{xi - y)DM{y - X2). 



(32) 



(33) 



We would like to stress that the representation (^) is completely equivalent to the initial representation (|l3|). The 
Green function G{xi,X2]X3^X4) can be considered as the Green function of the bilocal field A{xi,X2)- These fields 
are described by the nonlocal action S[A\. 



B. One-loop representation. 

The next problem is to give the standard particle interpretation to the action S[A\ in (^3|). For this aim this action 
should be represented in the form 



S[A] 



1 



{A+R-^A) + I,r,t[A] = 0(^3). 
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It means that we have to remove the term Unear in A and extract the quadratic term out of S[A]. Let us introduce 
the displacement 

^(2/1,2/2) 2/2) +^0(2/1,2/2) 



We get 



S[A] = -Eo - ^{A+A) - (A+Ao) - trln 



l + g(^A+^yD^]V 



(34) 



V = D 



M ■ 



1+5(A+x/a;;)£)m 



where the matrix multipUcation is implied. 

The constant term Eq is the vacuum energy in the lowest approximation: 

-^0 = -^{A+Ao) + trln [1 + 5 (^A+^/DZ) Dm 
and will be omitted in subsequent calculations. 



1. Linear term. 



The term linear in A should be equal to zero 

{A+Ao) + g tv[{A+./DZ)V] = 0, 



or 



Introducing the function 



we get the equation 



where 



Ao(a;i,a;2) = ^o(a;i - x^) 



a{x\ — 0:2) 
gsjDm{x\ - X2) 



a{xi - X2) = -g'^Dmixi - X2)'D{xi - X2) 



V = Dm- 



1 



1 



1 + o£>m' M2 + F + o(A;2) ■ 

Finally, we arrive at the equation of the Schwinger-Dyson type: 



~/,2^ 2 [ dp 



1 



)4 (m2 + (fc-p)2)(M2+p2 + a(p2))- 



(35) 



The integral in this equation contains the logarithmic ultraviolet divergence which can be removed by the renor- 
malization of the mass M . It means that we should put 

+ a(fc2) = + Urik"^, M^), 

where M^. is the "physical" mass of the constituent particle $ and 

a^(/e2,M2)=a(fc2)-a(-M2). 
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It is convinient to define the dimensionless function 



M2 



which satisfies the equation 



(m2 + {q- pf){M^ + p2 + M^w{p^)) 
1 



(36) 



{m? + {k-pf){M^ +p2 + M^w(p^))_ 
We have obtained the functional equation of the type 

w{k'^) = F{w,k\ 

Solution can be found by the fixed point method, i.e. we choose the initial "point" WQ{k'^) and calculate 

Wn+i{k'^) = F[wn,k% for n = 0,1,2,.... 

In the limit n ^ oo we get 

The main problem is to choose the zeroth approximation ti;o(fc^). We proceed in the following way. The propagator 
'D{k'^) after the mass renormalization should behave for k"^ — > like 



(37) 



Let us define the zeroth approximation 



l + c 



where the constant c is determined by the equation 

d 



dk 
M 



dp 1 



(27r)4 (m2 + {k- p)2) M2 + p2 



1 f da a{l — a) 



167rM2 TT J ^2Q,+ (l-a)2 



AZ-a(0. 



Here 



a(0 = - 

TT 



The equation 



-l-(l-r)ln^ + 



Z = 



arctan 



(38) 



l + c 1 + AZa(0 



gives 
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Thus, one can choose in the zeroth approximation 



In the first approximation (n=l) we have 



wi^-s(k^) ^ XZ — In -4 --T7 



where 



C(i) = 



= AZ(o)a(^) = c, 



Finally, we have in the zeroth and first approximations 



1 



2 for fc2 M^. 



^^^^ ^ M2 + fc2 + M2u;(i) (fc2) ^ Af2 + k 

One can check that 

Ao)(fc')«^(i)(fc') 

with 5% accuracy. 

The renormalized coupling constant gr can be defined by a standard way: 

5, = .gZ, A, - (41) 

because 

^^^'''^ " M2 + fc2 + MXfc2) m/+ fc2 f""^ 

In subsequent numerical calculations we use the zeroth approximation _D(o) ( pO| ) which gives quite acceptable 
qualitative semiquantative estimations. 

In this approximation the renormalized coupling constant is 

A, = A, (A, = = ^ '^^ (42) 

(l + v/l + 4Aa(e))2 

It is important that the renormalization considerably diminishes the initial coupling constant A. For example, in the 
case of the "deuteron", when <I>-particle is the "proton", 0-paricle is the "7r-meson" and A = 14.5 

Moreover, the renormalized coupling constant A^ is bounded for any A 

A. <^<^-15.01.... 

a(0 a(l) 
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S. From hilocal to local fields 



After removing the linear term we have 



S[A]=--{A+ A) -trim [l + g{A+ Dra)V 



hii(i + g) = in(i + Q)-g. 



The " trace" of " In" consists of terms 



dxjdxj+i !>{... - Xj)A'^{xj,Xj+i)^D„i{xj - Xj+i)'D{xj+i - ...) 



Let us proceed as follows. We introduce new variables 



- -¥1- 
Xj+i — Zj 2 I 



and the notation 

A{xj,Xj+i) ^ W{zj,yj), A{xj+i,Xj) = A+{xj,Xj+i) = W{zj,~yj). 
The term ( ^ ) can be written as 

J J dz,dy, D (... - z, ~ I) W{z,, -y,)^DMD {zj - \ 
dz, D{...~z,)V{z,,P,^)D{z,- ...) 



where 



V(z,P,) - / dy ^Dm{y)W{z,~y)e-'"-^' 



(43) 



(44) 



Px- 



dx — di 



Let the system of functions {C^q(2/)} with quantum numbers Q = (n/{/i}), where n, / and {/x} are radial, orbital 
and magnetic quantum numbers, be orthonormal, i.e., 

(45) 



(UqU*q,) ^ J y UQ{y)U*Q,{y) = 5qq> = 5nn-5u'5{p,}{,,'}, 

Y.UQ{y)U*Q{y')^5{y-y'), 
Q 

The function W can be represented by 

W{z, -y) = WQ{z)UQ{-y), W^{z) = Wq{z), 
Q 

W{p, -y) = WQ{p)UQ{-y), W^{p) = Wq{p), 
Q 



(46) 



Then we have 
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Q 
Q 



In this notation one obtains 

{A+A) = I dxi I dx2 A+ixi,x2)A{x2,xi) = {WW) 



Y^jdz Wq{z)Wq{z) = Y.j'^P Wq{p)Wq{p), 



tr Ini 



l+g{A+^/D^)V =trlni(l+5^(W)I?^) 



The basic representation for the Green functions under consideration becomes of the form 



G(xi,X2;x3,a;4)^GQ,Q,(zi,z2)- / \{DWq ■ Wq,{zi)Wq,{z2) ■ ^""^"^^ 



S[W] = ~^{WW) - trhii [1 + gr{WV)Vr] . 



C. Particle interpretation of the quadratic term. 



Let us extract the quadratic form from 

S[W]^-^{W[I- gMW) - trln2[l + 
ln2(l + s) = hi(l + s)-s + —. 

Here 

gfn = gMVVrVVr] 

and according to (|4^) and 

{Wg^.UW) 1 1 WQ{x)glXlQQ,{x - x')Wq[x') 

QQ' •' •' 

= Y.I dpW*Q{p)glnQQ,{p)WQ{p). 
Q 

The polarization operator g^Hqqi looks 

glnQQ.{x ~ x') = glir {VQVrVQ.Vr} 

= alVq (p,) Vr{x - x')Vq' (p,-) Vr{x' - x), 

^gl 1 1 dydy' UQ{y)P{x - x' ;y,y')U^,{y'), 



P{x; y, y') = v/A^y)^?. (a; - Vr [x + v/A^T), 

P(p;2/,y') = / dxe'^-P(:x-y,y'). 
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In the momentum space we get 



dk 



gtTlQQ,ip) = gt J VQ{k)Vr [k + P,. [k - VQ,{k). 

In our approximation we have 

dk VQ{k)VQ,{k) 



9r'^QQ'{p) = 9r 



(2^^)" (Af2 + (fc + f)^)(Af2 + (fc_£)^ 



(52) 



(53) 



The orthonormal system {Uq{x)} should be chosen so that the polarization operator IIqq'{p) should be diagonal 
in radial (n,n') and orbital {1,1') quantum numbers 



nQQ'(p) = Snn'Sw ■ Il["^^^^,^{p) 

The index structure of the diagonal polarization operator nj^j|^,j(p) looks like 

where the tensors ^\^}{^'}ip) contain combinations of the vectors p^p^'. 

The diagonal quadratic form of ( ^9[) gives the equation of motion for the field Wq{x) = ih}^^^ 



(54) 



d 



or 



Sqq' ~gliiQQ'{p)\ WQip)^Q. 
The requirement that this equation on the mass shell should be the Klein-Gordon equation gives the constraint 



(55) 



d 



on the mass shell. Thus, the function T^|^|'' (p) satisfies the equation 

\-glil^-'\p^)\w'CZ....^{p) = ^. 
The mass of the state with quantum numbers Q ~ (nl) is defined by the equation 

Let us write 

-1 + g'rfl^"'\p') - -Z(„,)(p2 + Aff„,)) + S("')(p2), 
E("')(p2)=g,^^(«^)(p^A^(^„,)) 

= 5.' W'^'Hp') - n("'n-Aff„o) - n'(„,)(-Mf„,))(p2 + Aff„,))' 



(56) 



(57) 



Z{nl) = gl -!!(„;) (-M, 



nl) / 



The constant Z(^ni) is positive. 

New field variables can be introduced as follows: 
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^q(-) - (58) 



The representation (Esh assumes of the form 



GQ,X^-y) = ^Qq,X^-V), (59) 



^Qo(a;-y)= / ¥'Qo(a;)'/'Qo(2/) 

where the action S looks hke 



-i(^p-V)-I^n^M (60) 



and the appropriate normaUzation should be chosen. 
The kinetic term is 



(^P- V) = + Ml + Eb] (p) (61) 



Q 



and the interaction term is 



Q 



X„t [(^] = tr ln2 [1 + (h^V)V\ , (62) 
(V^) = X! ^Q'^Q^Q^ = 



-n'Q(-A^^ 



The effective dimensionless coupling constants are defined as 



As a result, the final representation (|5^) can be interpreted as a partition function of the quantum field system of 
bosonic fields {0q} which have masses Mq and are described by the nonlocal action (|60|). 

We would like to stress that the resulting representation for the generating functional does not contain the initial 
coupling constant g. 

All calculations with the generating functional ( p9|) can be performed by perturbation expansions in coupling 
constsnts hq. We can trust these calculations if and only if the effective coupling constants (|6|) are small enough: 

x^A"^ « 1. 



D. The orthonormal system. 



The next step is to determine the orthonormal system ( [fSD . The problem is to find the spectrum and eigenfunctions 
of the operator P{p; y, y') in (|5l|), i.e. 

Pip)UQ ^ Eq{p^)Uq 



dy'V{p;y,y'}UQ{y\p)^EQ{p)UQ{y,p), Q = {n,l,{fi}). (64) 



20 



This equation can be represented in a standard form of the Bethe-Salpeter equation in the one-boson exchange 
approximation. Using the relation 

K+K^ ■ [ dx e'P^Vr (x - Vr (x + = 5{y - y') 



with 



and introducing the functions 



we get the standard form of the Bethe-Salpeter equation (see, for example, 



dy 



*Q(y,p) =5r^m(y)*Q(y,p), 



where the spectrum is defined by the equation 



(65) 



(66) 



Thus the diagonalization of the operator P{p;y,y') is equivalent to the solution of the Bcthe-Salpeter equation in 
one-boson exchange approximation. 

It should be noted that in this formulation the functions 5'Q(y,p) depend on the direction of the vector p — (pojP) 
in the space R'', the standard choice is Pni = (iM„;, 0). 

The main problem to use the Bethe-Salpeter basis is that the Bethe-Salpeter equation can only be solved by 
numerical methods. Even the solution obtained by Wick and Cutkosky is reduced to the differential equation 
which should be numerically computed. Our aim is to continue analytic calculations as long as possible in order 
to get a visible general picture of arising bound states in the system under consideration. Therefore we choose a 
more practical way, namely, we use an orthonormal basis that is simple enough from an analytic point of view and is 
directly connected with the problem under consideration. In this case the operators g^T^QQ' are not diagonal so that 
we should diagonalize them. The idea consists in finding an effective basis for diagonalization of g^Hqq' such that 
its lowest function would provide a good qualitative description for the eigenvelues and the next two or three 

functions only give a good quantitive description for those eigenvelues. 

This effective basis {Uq{x)} can be constructed using the boson Green function Da{u) with a parameter a as a 
weight function inducing uniquely the system of orthonormal polynomials in the space R^. Thus, the full orthonormal 
system of functions (HS) can be chosen in the form 



UQ{x,a) = Da{x)aPQ{ax). 

Here Pq{u) are real polynomials, satisfying 

Pq{-u)^{-1)'Pq{u), so that U*q{x) = Uq{-x). 

The parameter a that enters into the orthonormal system is fixed by a variation condition formulated below. 
The construction of this basis is presented in the Appendix A. 



(67) 



E. The diagonalization procedure. 

Let us demonstrate our procedure of diagonalization of the polarization operators (|5^) for states Q — (nO) and 
Q' = (n'O). In the momentum space this polarization operator looks like 



dk V(no) {k\m, a)V(n'o) {k\m, a) 

(2^)' ' (m2 + {k+ f )') (Af2 + (fc - f )' 



Srii(«n')(P ) - 9r TTTu ' 7 TtTI ^ 72\ ^^^^ 
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where the vertex defined by (E^) looks in this case like 



V(nO)ip\m,a) ^ dx ?7(„o)(a;)\/-D™(a:)e'^ 



ixp 



(69) 



= dx Dm{x)Da{x)aP(^nO){ax)e''-'P . 



One can develop 



yjD„,,{x)DmAx) - - ^™'(-^) = 



with m = 2^4^ and A = JUi^ and 



We shall use the lowest approximation 



n=0 



A 



A2 
m 



dk 



(27r)4 (to2 + /c2)l+r, 



^Dra{x)Da{x) « (a;), 
the accuracy of which is quite acceptable for our consideration. Then, the vertex acquires the form 



(^) +p2 



(70) 



(71) 



The explicit form of the polynomials P(„o)('"^) is given in the Appendix A. We get 

1 



R(n){t) 



dt R(^n){t)R{n'){t) 



'1 



Abt{l - t) 
(1 - htf 



1 



,1 -t 



rn, a 



Now we formulate the variational principle which defines the parameter a. The operator 11(00) is the largest 
eigenvalue of the operator matrix !!(„„/■), therefore the parameter a can be defined by the variation requirement 



max gl 11(00) (-A/.^^,) 



(72) 



— max — / dt 

v 27r 



i±27 

2 



t+l-t 



Abt{l - t) 



{1-htf 



- 1 



where the notion 



(0) 



2Mr 



m a 



are used. Thus, the parameter a = a{M^Q-^, m) is a function of m and A/(o). Quite a good approximation is as follows 

(73) 



r/ = 77(e,6)«2-^6+^e. 



Let us show that this orthonormal functions with the parameter a ( fz^ ) gives quite good approximation for the 
eigenvalues of the matrix !!(„„/) (p^^ por this aim we calculate the matrix 



r 



N = |n(„„') (6,^,77), {n,n' = 0, ...,iV)| 
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and their eigenvalues 



f^ = diag{i;r,...,i;r}. 

Then we have to compare Ej^^ for fixed j and difi'erent N. 

The numerical results are given in Table I. The first case for ^ = .5, b = .25, r] = 2.451 and the second case for 
^ = .2, b = .9, r/ = 1.22. One can see that for the lowest eigenvalue practically the first lowest eigenfunction can be 
used, i.e. our choice of the orthonormal system gives quite a good accuracy. 



Table I. Diagonalization of the matrix Vj\f. 



N 




El 


E2 


E3 





.04165 








1 


.04166 


.009941 






2 


.04173 


.010279 


.002755 




3 


.04175 


.010368 


.003295 


.0007482 


4 


.04175 


.010402 


.003546 


.0010336 





.1239 








1 


.1262 


.03564 






2 


.1262 


.03616 


.01162 




3 


.1263 


.03645 


.01298 


.003789 


4 


.1263 


.03655 


.01373 


.004710 



F. Applicability of S-, BS- and BNC- methods. 

Three parameters ^, A,- and b are not independent. The standard formulation of the problem is the following: the 
parameters ^ and are given and the mass of a bound state b has to be found. We reformulate this problem: what 
is the region of changing b for the fixed parameter ^, if the effective coupling constant Ae// is smaller then one? 

In this section we give qualitative curves 

AM M(o)(0 , 

2M 2Mr 

which restrict the admissible regions of ^ and b when the effective coupling constant in S-, BS- and BNC-methods 
is smaller then 1. These curves are presented in Fig 1. and show the applicability of S-, BS- and BNC-methods to 
study the bound state problem. 
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FIG. 1. Applicability of S-, BS- and BNC- methods. 

The region of applicability of the nonrelativistic Schrodinger equation is a small vicinity near point (0, 0) in 
Fig 1. 

The range of existence of bound states (under solid line in Fig 1) is defined by the equation 

l = A„a.(C)n(o)(C,6) 

where 



1 



n(o)(^,&) 



dt 

max 
v J ^^t: 





t + l-t 



The region of applicability of the Bethe-Salpeter method (under dashed line in Fig 1) is defined by A < 1 and 

A. <Ai(0 = A,(l,0 ^ 



(1 + 



where b — 6(f) is defined by the equation 



i = Ai(e)n(o)(e,6(0)- 



The rigion of applicability of the BNC-mcthod (under dotted line in Fig 1) is defined by the inequality 



where the effective coupling constant is defined by ([63 



1 
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G. Small binding energies. 



In this section we show that the BNC-method gives the correct nonrelativistic hmit in the case of the Coulomb 
potential. For the Yukawa potential we show how the screening length depends on the coupling constant and the 
patanieter £, = If the binding energy is small then 



e = 1-6: 

and one can get for (^4) when e is small 



\2M 



2M - Mh AM 



M 



M 



< 1. 



1 

2^ 



dt 



{^) t + l-t 



4bt{l - t) 



(1 - btf 



V(^) (1-^^) 



eu 
2 



4w(l - eu) 



e3/27r 



duy/u 



e(l + u)2 
477^ 



l + « (e + ?7)(e + ^ + 2V^)2- 



The Coulomb potential H — ^ in our approach is defined by the conditions m = or ^ = and a ^ 1. We 



get 



and from 



it follows that 



11(0) (0, = max 



4?7 



1 



n (ry + 27^)2 2yi 



l = A,n(o)(0,e), or 1 



2x/i' 



{Xr — a) 



^ , , Mr o? ^, 

AAf = = Mem 

2 2 2 

according with the lowest energy of the nonrelativistic Schrodinger equation for the Coulomb potential. 

2 

For the Yukawa potential H ~ ^g-mr -j-j^g critical screening length, defined by the condition e = for the 

lowest state, can be computed. For small ^ this length is defined by the equation 



1 — a max 



4j^2 



16 a 



and 



32 Mr 

rUc = 27"~2~ ^ 1-1851... • aMcm- 
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32 



1.1851. 



aMcn, 27 
The numerical calculations |l^ give 

rUc = 1.1906... • aMcm- 

In the case < ^ < 1 the relation between ~ ^ and ttc is given by 

Uc ■ T{^c) = 1 and = 2f,r(Ce) 

where 

1 / 



T(ec) = max / 



1 



At 



- 1 



The results are given in Table II. One can see that for > -1 the relativistic corrections are more then 5%. 

Table II. The critical screening length for the Yukawa potential. 



^ Mr 






0.0001 


1.1851 


0.00016875 


0.001 


1.1846 


0.0016888 


0.01 


1.1795 


0.01695 


0.1 


1.1324 


0.1766 


0.2 


1.0861 


0.368 


0.4 


1.0079 


0.793 


0.6 


0.9436 


1.3735 


0.8 


0.8893 


1.799 


1.0 


0.8426 


2.373 



V. APPENDIX A. ORTHONORMAL SYSTEM 

The orthonormal system {Uq{x)} is defined by 



Uq{x) = i'^ Da{x)aPQ{ax 
a" 1 



Da{x) 



(27r)2 ax 



Ki{ax), X = V x' 



Q = (n/{^}), {^} = 

d'^x Uq{x)Uqi{x) = I d'^x Da{x) aPqiax) aPQ'{ax) = Sq 



d'^x — j drr^ j dn, J dn = J '^'t' J '^^ ^^"^ 9 J dj sin^ 7 
000 
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The polynomials are defined by 



X u 

u = ax, n - 



The orthonormality condition looks as 

dx Da{x)aPQ{ax)aPQ'{ax) 



= . " / dxx^ ■ —Ki{ax) ( dn PQ{ax)PQi{ax) 
{My J ax J 



oo 

^ y duu^ -f^i(w) j dn Pq{u)Pqi{u) 





(27r 



(2§^ ■ ■ 2'(r+l) '^"' ■ ^W^'"'} " ^^'5'' " 22x/r+l, 

OO 

^Jds S^+^'K^{s)P(^i){s^)P(^.i){s^)=5nn' 



A. Angular polynomials T'^^-^{n) 



Polynomials T^^^j (n) satisfy 



• is symmetric for /Xj ^ iij, 

• Tjil.^,{n) = iP(l|2..i)T^'-^),(n) - ^^P{l2\Z...l)5,,^,Tg-.^l{n). 
The first four polynomials are 



^123 W = nin2ns - ^(ni^23 + ^2^31 + n3^i2), 
^i234(") = nin2n3n4 

--{nin2dsi + nin3(524 + nin4d23 + n2n3^i4 + n2n4e/tai3 + 713714^12) 

o 

+ 7^(^12^34 + <5l3^24 + ^14<523) 
4o 



where the condensed notation Uj = n^^ and 6ij = 5^^. are introduced. 
The normalization of the angular polynomials is defined by 
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{m} 



dt \J\-t'^ ■pi{t)pi,{t) = / d7 sin^7 (cos 7)/?;/ (cos 7) = 



where (t) are the Gegenbauer polynomials. In addition one can get 

dn (T('),(ni)T('),(n)) • iT['hi^)T^'hi^,)) 



2'(? + l) 



22'(; + i) 



B. Radial functions 



The radial functions are defined in a standard way (see, for example, p9|). Let us define 



c{n) = ^ I du u^+^"Ki{u) ^ 4"(n + l)!n!. 



CniD = 



c{l) c{l + l) ... c{l + n) 
c{l + l) c{l + 2) ... cil + n + l) 

c{l + n) c{l + n + l) ... c{l + 2n) 



The radial polinomials read: 



Pinl){ 



V£'n-l(0^n(0 
c(/) C(/ + 1) C{1 + 2) 

c{l + l) c(/ + 2) cll + 3) 
c{l + n-l) c{l + n) c{l + n+l) 



1 



where Z?„(Z) = det Cn{l). 

The orthonormality condition is as follows: 



c{l + n) 
c{l + n+l) 

c{l + 2n- 1) 
(u2)» 



1 



du V^^"^^ Ki{u)Pi^ni){u'^)P(n'l){u'^) = 5nn'- 



The first four polynoms are 



aP{m) = a, 
aP(w){ax) = 



aP{20) {ax) 
aP{30) {ax) 



-1 
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a 



ba^x"^ {a'^x^)^ 
3 — + ^ ^ 



96 



V4303 



-34 



370^x2 53(0^0;^)^ 13{a^x^)^ 



192 



9216 
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C. Vertex functions 



The vertex function 

Vqik) =t' Jdx ^Dm{x)Da{x) aPQ{ax)e-''''' 

in the approximation 



^/D^{x)Da{x) « Dr^ (x) 

looks as 



VQ{k) ^i^Jdx D^{x)aPQ{ax)e-'^^ = i^Pg (za^^ ■ ^ 



The following equality is useful: 



\dk^) ' + k^ + k^ 



j=0 



4"s!(s + l)!(s + j)! (fc^V 

+ l)l{s - f)l ' (M2 + A:2)«+i+i 
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